
Figure 2: 



29 



00 



0.93 



0.0 



log* 



4.0 



Figure 3: 




Figure 4: 



30 




Figure 5: 



31 




Figure 7: 



32 



0.0 



logt 4.0 



Figure 8: 




Figure 9: 



33 




Figure 10: 



34 




Figure 12: 



35 



0.8 



0.0 log* 4.0 



Figure 13: 




Figure 14: 



36 




Figure 15: 



37 



Finite viscoelasticity of filled rubbers: the effects of 
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Abstract 

Constitutive equations are derived for the viscoelastic behavior of filled elas- 
tomers at isothermal loading with finite strains. A particle-reinforced rubber is 
thought of as a composite where regions with low concentrations of junctions be- 
tween chains are randomly distributed in the bulk material. The onset of these in- 
clusions is associated with the inhomogeneity in spatial distribution of a cross-linker 
during the mixing process. With reference to the theory of transient networks, the 
time-dependent response of an elastomer is modelled as thermally activated pro- 
cesses of breakage and reformation of chains in domains with low concentrations of 
junctions, whereas junctions in the bulk medium are treated as permanent. Stress- 
strain relations are developed by using the laws of thermodynamics. Adjustable 
parameters in the constitutive equations are found by fitting experimental data in 
tensile relaxation tests for several grades of unfilled and carbon black filled natural 
rubber. It is demonstrated that (i) the average relaxation time noticeably grows 
with the elongation ratio, which is explained by mechanically-induced crystalliza- 
tion of strands, and (ii) the relaxation spectrum of a filled elastomer is not affected 
by mechanical pre-loading and thermal recovery at elevated temperatures. 
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1 Introduction 



This paper is concerned with the viscoelastic behavior of unfilled and particle-reinforced 
rubbers at isothermal loading with finite strains. The time-dependent response of filled 
elastomers has attracted substantial attention in the past decade jlL |2], |5], [13], [M], H| [T^, 
^ |26|, [27], [32], |33]. [34]. [3^. [36], |3!| , which is explained by numerous industrial 
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applications of these materials (vehicle tires, shock absorbers, earthquake bearings, seals, 
flexible joints, solid propellants, etc.). 

Observations in uniaxial and biaxial tests evidence that rubbers reinforced with carbon 
black or silica, as well as solid propellants filled with hard particles reveal strong time- 
dependent response 0, 
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3|,|39|. 



m m m m 

A particle-reinforced rubber is conventionally thought of as a network of macro- 
molecules bridged by junctions (chemical and physical cross- links, entanglements and 
filler particles). The classical theory of rubber elasticity |3l| treats all junctions as per- 
manent, which implies that chains cannot detach from the junctions. To predict the 



time-dependent response of rubbery polymers, Green and Tobolsky W5\ presumed that 
chains can slip from the junctions as they are thermally agitated. The concept of tem- 
porary junctions was developed in |2^, [37], f4(| and was applied to the description of the 
viscoelastic behavior of rubbers in ||, fTl] . 

We treat an elastomer as a strongly inhomogeneous material and assume that the 
characteristic length of an inhomogeneity substantially exceeds the radius of gyration for 
a macromolecule, on the one hand, and it is noticeably less than a size of a specimen, on 
the other. 

To describe the viscoelastic response of a rubber (unfilled or particle- reinforced) , we 
model it as a composite medium, where regions with low concentrations of junctions 
(RLCJ) are randomly distributed in the bulk material. These domains arise during the 
preparation of rubber samples (at the stages of mixing and vulcanization) because of 
the spatial inhomogeneity in the distribution of a cross-linker. We treat junctions in the 
host matrix as permanent and associate the time-dependent response of elastomers with 
breakage and reformation of chains in inclusions with low concentrations of junctions. 
The qualitative difference in the behavior of chains in RLCJs and in the bulk medium 
is explained by the fact that surrounding macromolecules in the bulk material prevent 
slippage of chains from the junctions, whereas they do not affect the slippage process in 
RLCJs, where interactions between chains are noticeably weaker. 

According to the theory of temporary networks |37f] , a strand whose ends are linked to 
separate neighboring junctions is treated as an active one. Snapping of an end of a strand 
from a junction is thought of as its breakage (transition from the active state of the strand 
to its dangling state). When the free end of a dangling strand captures a junction (which 
reflects creation of a new physical cross-link or entanglement), a new active strand merges 
with the network. Breakage and reformation of active strands are treated as thermally 
activated processes: attachment and detachment events occur at random times as the 
strands are agitated by thermal fluctuations. 

The objective of this paper is three-fold: 

1. to derive constitutive equations for the time-dependent behavior of particle- reinforc- 
ed rubbers using the concept of temporary networks, 
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2. to carry out uniaxial tensile relaxation tests on virgin specimens of unfilled and 
particle-reinforced natural rubber, as well as on the same specimens after mechanical 
pre-loading and recovery at an elevated temperature, 

3. to apply the stress-strain relations to the analysis of experimental data. 

Our purpose is to evaluate the effects of pre-loading and thermal recovery on the average 
rate of stress relaxation and on the relaxation spectrum. 

The study is organized as follows. Breakage and reformation of active strands in 
the regions with low concentrations of junctions are discussed in Section 2. Governing 
equations for the kinetics of rearrangement are introduced in Section 3. The strain energy 
density is determined in Section 4 for the bulk material, in Section 5 for RLCJs and in 
Section 6 for the composite medium. Stress-strain relations at finite strains are derived 
in Section 7 by using the laws of thermodynamics. Uniaxial extension of a specimen is 
analyzed in Section 8. Section 9 deals with the description of tensile relaxation tests. 
Adjustable parameters in the constitutive equations are found in Section 10 by fitting 
observations. A brief discussion of results is presented in Section 11. Some concluding 
remarks are formulated in Section 12. 



2 Breakage and reformation of active strands 

Breakage and reformation of strands in RLCJs are modelled as thermally activated pro- 
cesses. Slippage of an end of an active strand from a junction and reformation of a 
dangling strand occur at random times when they are thermally agitated. With reference 
to the theory of thermally activated processes [12[], the rate of breakage for an active 
strand in a stress-free medium is given by 

r = r - exp (-^r)' « 

where Co > is the energy for breakage, ks is Boltzmann's constant, and the pre-factor T* 
weakly depends on temperature. Equation (1) implies that for any u), the rate of breakage 
vanishes at low temperatures (T — > 0) and it approaches at elevated temperatures 
(T-oo). 

The dimensionless energy for breakage reads 

UJ 



knT° 



where T° is some reference temperature. Disregarding the effect of small increments of 
temperature, AT = T — T°, on the rate of breakage, T, we find the rate of transition from 
the active state of a strand to its dangling state 

T = T^exp(-u;). 

In what follows, this equation is assumed to be satisfied for an arbitrary loading program, 
provided that the coefficient depends on current strains and some internal variables 
that characterize 
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1. damage of the secondary network and isolated clusters of filler, 

2. changes in the structure of macro molecules driven by their alignment and mechanic- 
ally-induced crystallization. 

This implies that for a time- dependent loading, the rate of breakage for active strands 
reads 

r(t,w) = r»(t)exp(-w), (2) 

where the function T*(t) is found by fitting observations. 

Denote by p(ou) the distribution function for energies of breakage for active strands 
in RLCJs. The function p{uj) is assumed to be robust with respect to the influence of 
thermo-mechanical factors. This means that p is a function which does not depend on 
temperature (for small increments, AT) and on the history of loading. To fit experimental 
data, we apply the quasi-Gaussian distribution 

(3) 

where Q and S are adjustable parameters, and the constant po is found from the condition 

/ p(u)dw = 1. (4) 
Jo 

As it will be demonstrated later, the ratio 




that characterizes the width of the quasi-Gaussian distribution (3), is relatively small 
(about 0.5), which implies that Q may be associated with the average energy for breakage 
of an active strand, and £ determines the standard deviation for energies of breakage. 

The function r*(t) in Eq. (2) and the parameters Q and X in Eq. (3) entirely describe 
slippage of active strands from temporary junctions in the regions with low concentrations 
of cross-links. 



p(u) = p exp 



2E 2 



3 Kinetic equations for rearrangement of strands 

The kinetics of breakage and reformation of active strands is determined by the function 
X(t, r, oj), which equals the current (at time t) number (per unit mass) of active strands in 
RLCJs with the energy for breakage u that have last been bridged to the network before 
an instant r G [0, t] (the initial time t = corresponds to the instant when external loads 
are applied to the specimen). 

The quantity X(t,t,u) equals the number of active strands (per unit mass) with the 
energy for breakage oj at time t. In particular, X(0,0,cj) is the initial number of active 
strands with the energy for breakage u, 

X(0,0,w) = Nu>(u), (5) 
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where is the average number of active strands per unit mass of RLCJs. The amount 

dX 



dr 

t=T 



is the number (per unit mass) of active strands with the energy for breakage u that have 
been linked to the network within the interval [r, r + dr]. The quantity 

dX f 

— (t,T,Lu) dr 

determines the number of these strands that have not been broken during the interval 
[t, t] . The amount 

-—(t,0,u)dt 

is the number of active strands (per unit mass) that detach from the network (for the 
first time) within the interval [t,t + dt], and the quantity 

9 2 X . . , , 
-{t, r, uj) dtdr 



dtdr 



is the number of strands (per unit mass) that have last been linked to the network within 
the interval [r, r + dr] and slip from temporary junctions (for the first time after merging) 
during the interval [t, t + dt]. 

The evolution for a transient network is determined by the relative rate of breakage for 
active strands, T(t, uj), and by the rate of merging with the network for dangling strands, 
j(t, uj). The assumption that Y and 7 depend on the energy for breakage, u, distinguishes 
the present model from previous ones. The dependence of the rates of breakage and 
reformation on time, t, reflects the effects of strains and mechanically-induced changes in 
the internal structure of chains on these parameters under a time-dependent loading. 

The rate of breakage, T, equals the ratio of the number of active strands broken per 
unit time to the total number of active strands. It follows from this definition that for 
strands connected with the network before loading and remaining active at time t, 

1 f)X 

Applying the definition of the rate of breakage to strands merged with the network at 
time t > and remaining active at time t > r, we find that 



r(t,u) = - 



dX, 



-1 -f)2 



d 2 X 
dtdr 



(t,r,u). (7) 



It is assumed in Eqs. (6) and (7) that the rate of breakage, T, is independent of the 
instant, r, when an active strand merges with the network. 

The rate of reformation, 7, equals the number of dangling strands (per unit mass) 
bridged to the network per unit time 



, \ dX, 



(8) 

T=t 
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Integration of differential equations (6) and (7) with initial conditions (5) and (8) implies 
that 



X(t,0,u) 
dX 



dr 



(t,T,U>) 



N L p(u) exp 
7(r, a;) exp 



— / r(s,ou)ds 
Jo 

— J T(s,uj)ds 



(9) 



To exclude the rate of reformation, 7, from Eq. (9), we suppose that the concentrations 
of active and dangling strands are independent of time. This means that the number of 
strands broken per unit time coincides with the number of strands merging with the net- 
work within the same interval. The number of active strands detached from the network 
during the interval [t, t + dt] equals the sum of the number of initial strands (not broken 
until time t) that slip from temporary junctions 



dx 

~dt 



(t,0,cu) dt 



and the number of strands linked with the network within the interval [r, r + dr] (for the 
last time before instant t) for various instants r and broken within the interval [t, t + dt] 



d 2 X 
dtdr 



(t, r, a;) dtdr. 



This number coincides with the number of strands, "f(t,iv) dt, merged with the network 
within the interval [t,t + dt], which results in the balance law 



7(f,w) 



dX 



(t,0,Uj)- I 

Jo 



* d 2 X 



dt Jo dtdr 

Substitution of Eq. (9) into this equality implies that 



(t, t, uj)dr. 



-f(t,u) = T(t,uj)^N L p(uj) exp 

The solution of the linear integral equation (10) reads 

7 (t,w) = N L p(cu)T(t,u) 



J T(s,uj)ds + J 7(7-, cj) exp — J T(s,uj)ds rfrj. (10) 



(11) 



Equations (9) and (11) determine the current distribution of active strands in the regions 
with low concentrations of junctions. 



4 Mechanical energy of the bulk material 

The bulk medium is treated as a permanent network of long chains bridged by junctions. 
We adopt the affinity hypothesis which disregards thermal oscillations of junctions and 
presumes that the deformation gradient for the motion of junctions at the micro-level 
coincides with the deformation gradient for the motion of appropriate points of a sample 
at the macro-level. 
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The bulk material is modelled as an isotropic incompressible medium, whose strain 
energy density, Wb, depends on the current temperature T (the strong effect of temper- 
ature reflects the entropic nature of the mechanical energy for elastomers), a vector of 
internal variables a (which accounts for changes in the structure of a particle-reinforced 
rubber driven by rupture of filler clusters, stress-induced alignment of macromolecules, 
and their crystallization) and the first two principal invariants Ik (k — 1,2) of the Cauchy 
deformation tensor 

C°(f) = V°R(t) • [v°R(t)] , 

where V°R(£) is the deformation gradient at time t, V° is the gradient operator in the 
reference state, the dot stands for inner product, and T denotes transpose. The strain 
energy density reads 

W B (t) = W B (T(t), a (t),h(C°(t)),I 2 (C°(t))), (12) 

where the function Wb obeys the condition 

W B (T, a ,h,I 2 ) =0. (13) 

il=d,i2=J 

Equation (13) means that the mechanical energy vanishes in the stress- free state. In 
this equation we disregard thermally induced changes in the reference state of rubber 
compared to the strains under consideration. 



5 Mechanical energy of RLCJs 

Regions with low concentration of junctions are assumed to be randomly distributed in 
the bulk material as inclusions with various sizes and shapes. Any RLCJ is modelled 
as an isotropic and incompressible transient network, where breakage and reformation of 
active strands are driven by thermal fluctuations. The strain energy of an active strand 
in regions with low concentration of junctions, w^, is assumed to depend on temperature, 
T, the vector of internal variables, a, and the first two principal invariants of the Cauchy 
deformation tensor. The conventional hypothesis that stress in a dangling strand totally 
relaxes before this strand captures a new junction implies that the stress-free state of a 
strand merging with the network at time r > coincides with the deformed state of the 
network at that instant. This means that for a strand not broken within the interval [0, t], 

w L (t,0) = w L (T(t),a(t),I 1 (C°(t)),I 2 (C°(t))) } 

and for a strand that has last been reformed at time r e [0, t], 

w L (t, r) = w L (T(t), a(f), h(C(t, r)), I 2 (C(t, r))) , 

where C(t, r) is the relative Cauchy deformation tensor for transition from the deformed 
state of a specimen at instant r (when the strand has last been merged to the network) 
to its deformed state at time t > r. The tensor C(t, r) is expressed in terms of the 
deformation gradient V(r)R(t) by the formula 

C(t,r) = V(r)R(t)-[v(r)R(t)] T , 
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where V(r) is the gradient operator in the deformed state of the medium at time r. 

Summing mechanical energies for active strands with various energies for breakage, 
u, merged with the network at various instants, r G [0,t], and neglecting the energy 
of interaction between strands, we arrive at the strain energy density per unit mass of 
RLCJs, 

roo r rt QX 

W h (t)= X(t,0,u)w h (t,0)+ —(t,T,u)w L (t,r)dT du. 
Jo L Jo or 

Substitution of expressions (9) and (11) into this equality implies that 

f OO f ft 

W L (t) = J p(u)du^W L (T(t),I 1 (C (t)),I 2 (C (t)))exp[-j o T(s,uj)ds 



+ jT* r(r, u)W L (T(t), h(C(t, r)), J 2 (C(t, r 



x exp 



— / T(s,uj)ds 



dr 



where the function W L (T : I 1: 1 2 ) = N L w L (T, Ii, I 2 ) satisfies condition 



W- L (7\a,Ji,J 2 ) 



Ji=3,/ 2 =3 







(14) 



(15) 



similar to Eq. (13). Formula (14) determines the strain energy density (per unit mass) of 
inclusions with low concentrations of junctions. 



6 Strain energy of a filled rubber 

Denote by /il the mass fraction of regions with low concentrations of junctions. Taking 
into account that the strain energy density of a composite material, W, is an extensive 
thermodynamic variable, we find that 



W = (1 - fi L )W B + fi L W L . 
Substitution of expressions (12) and (14) into Eq. (16) results in the formula 

wit) = (l-^w^m^^j^it^hic^t))) 

+/i L J™ p(u)du{w L (T(t), a(0, h(C°(t)), h(C°(t))) exp 
+ J* r(r, u)W L (T(t), a(t),h(C(t, r)), J 2 (C(t, r))) 



(16) 



r(s, u)ds 



x exp 



T(s, u)ds 



dry 



(17) 



It follows from Eqs. (15) and (17) that 



^(t) = ^)^W+A(t)-^(t)+f/rW^(c o (t))+f/ 2 °(t)^(c°(t))+f/(t)-r(t), (is) 
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where 
J(t) 



(l-^)^(T(t),a(t),h(C°(t)),I 2 (C°(t)) 

roc ( dWi / \ r /"* 

+// L ^ pM^[^(T(t),a(t),/ 1 (C°(t)),/ 2 (C (t)))exp[-y o r(s,u;)ds 



dWr 



dr 



A(t) 



+ I r(r,c)^(T(t),a(t),/ 1 (C(t,r)),/ 2 (C(t,r)))exp[- jf_ I^ds 
(1 " ML) ^( T(t) ' a(t) ' /l(C0(t)) ' /2(C ° (t)) ) 

+/i L ^ pM^[^(T(t),a(t),/ 1 (C (t)),/ 2 (C (t)))exp[-y o r(s,u;)ds 



U(t) 



+ / r(r,o; 



dW L 
da 



T(t),a(t),/ 1 (C(t,r)),/ 2 (C(t,r)) exp 



T(s, ui)ds 



dr 



{i- ^)^(T{t)Mt),h{c o {t)lh{c°{t)) 



dh 



/i L / o 00 ^(T(t),a(t),/ 1 (C (t)),/ 2 (C°(t)))exp[-| o r( SlW )rf S ]pH^, 



p(u)du / r(r, u;) exp 



r(s, u)ds 



dr 



x 



dW L 



flr 

(T(t), a(i), J^Cft r)), h(C(t, r))) ^(C(i, r)) 



Of 



dh 



dh 



(T(t), a(t), I^t, r)), I 2 (C(i, r))) -^(C(t, r)) 



y(t) 



/i L y o r(t,^M|iy L (r(t),a(t),/ 1 (c o (t)),/ 2 (c°(t)))exp[- ^ r(s,u)ds 

+ f r(r, o;)W L (T(*), a(t), ^(Cft r)), J 2 (C(t, r)) 



x exp 



T(s, u;)<is 



\<ko. 



(19) 



The principal invariants of the Cauchy deformation tensor, C(t, r), coincide with the 
principal invariants of the Finger tensor for transition from the deformed state at time r 
to the deformed state at time t, 



F(f,r)= V(r)R(t) • V(r)R(t). 



This assertion implies that 



^(C(t l r)) = ^(F(t l r)) = §(F(t l r)): 



<9F ^ T 



(20) 



where the colon stands for convolution of tensors. The derivatives of the principal invari- 



ants, K, read 



dh 
dF 



dh 
dF [ 

9 



h(F)I — F 1 



where I is the unit tensor. The derivative of the tensor F with respect to time is given 

by § 

^(t,r)= [v(t)V(t)] T -F(t,r) + F(t,r)-V(t)V(t), (21) 

where V is the velocity vector. It follows from Eqs. (20) and (21) that 

dh 

dt 
dh 



;C(t,r)) = 2F(t,r):D(t), 



dt 



[C(t,r)) 



h(C(t,r))F(t,r)-F 2 (t,r)\ : D(t), (22) 
is the rate-of-strain tensor. Substitution of 



where D(t) = §[V(t)V(t) + (V(i)V(t) 
expressions (22) into formula (19) for the functional U implies that 

U(t) =2/i L A L (t) :D(t). 

Here 



(23) 



t r 



p{u)du I L1 (t,r)F(t,r) +0 L2 (t,r)F i (t,r) 
rt 

dr, 



o 



xT(r, u)exp — J T(s,uj)ds 
and the functions 4>Lk(t,r) read 



9/i 



0L2(t,T) 



+Ji(C(*, r))^(r(*), a(f), Ji(C(t, r)), I a (C(t, r)) 
9/ 2 



(T(t),a(t),/ 1 (C(t,r)),/ 2 (C(t,r 



By analogy with Eq. (23), we find that 

dh 
dt 



XX°(*)#(C°(t)) = 2{(1-/, l )AbW 



fc=i 



+H h A° h {t) J™ exp[- J*T{s,uj)ds\p{uj)duj} : D(t), (24) 



where 



A;(t) = 0° (t)F°(t) + 0° (t) F°(t) G9 = B,L) 



F°(i) = V°R(t) • V°R(t) is the Finger tensor for transition from the reference state of 
the medium to its deformed state at time t, and the functions 4>°g k {t) are given by 



tin® 



dh 



T(t), a(t), h(C°(t)), I 2 (C°(t)) 
dWi 



<t>°82{t) 



+h(C°(t))^(T(t),a(t),h(C°(t)),h(V 
9Wp {T(t),a(t),h(C°(t)),I 2 (C°(t))). 
10 



dh 



It follows from Eqs. (18), (23) and (24) that 

™{t) = J(t)^(t) + A(f) • ^(f) + 2H(f) : D(t) - Y(t), (25) 

where 

H(t) = (1 - /hJAbW + A^(t) exp [- J T{s, u)ds]p{u)du + A L (t)}. 

7 Constitutive equations 

We confine ourselves to quasi-static loadings with finite strains, when mechanically in- 
duced changes in temperature, AT, are rather weak. This implies that the effect of tem- 
perature on the specific heat, c, may be neglected. The following expression is adopted 
for the free energy (per unit mass) of a filled elastomer: 

V = + ( c - S°)(T -T°) -cT\n^ + W, (26) 

where \1/ and S° are the free energy and the entropy per unit mass in the stress-free state 
at the reference temperature T°, and the second and third terms on the right-hand side 
of Eq. (26) describe the free energy of thermal motion of chains. 

For an incompressible medium, the Clausius-Duhem inequality reads || 

dQ dT d-% 1/ . 1 \ 

where p is mass density, X' is the deviatoric component of the Cauchy stress tensor S, q 
is the heat flux vector, S is the entropy per unit mass, and Q is the entropy production. 
Substituting Eq. (26) into Eq. (27) and using Eq. (25), we find that 



P 



+A/ ft +Y -^ T> -°- (28 > 

Because Eq. (28) is to be satisfied for an arbitrary loading program, the expressions in 
brackets vanish. This implies the formula for the entropy per unit mass 

S(t) = S° + c\n^--J(t) 

and the constitutive equation 

E(t) = -P(f)I + 2p{(l-/i L )A£(f) 

+/i L [AS? (f) J™ exp(- J* r(s, u)ds)p{u)dw + A L (t)] }, (29) 

where P is pressure. 
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We adopt the Fourier law for the heat flux vector q, 



q = - K VT 

where k > is thermal diffusivity. This implies that the last term on the right-hand 
side of Eq. (28) in non-negative. It follows from Eq. (19) that the functional Y is non- 
negative as well. To ensure that the Clausius-Duhem inequality (28) is satisfied, it suffices 
to assume the derivative of a to be proportional to the vector A, 

da 

it = ' 

where rj > is a function to be determined. 

To distinguish explicitly rearrangement of temporary networks in RLCJs and rupture 
of filler aggregates, we suppose that i] is a function of the vector a, the rate-of-strain 
tensor D, and the Cauchy deformation tensor C°, t] = r](a, D, C°), which vanishes for 
any time-independent straining of a specimen, 



77(a,D,C c 



= 0. 

D=0 



In the sequel, we concentrate on the study of tensile relaxation tests and do not dwell on 
the kinetic equations for the vector of internal variables. 



8 Uniaxial tension of a specimen 

We apply constitutive equation (29) to analyze stresses in a bar under tension. Points 
of the bar refer to Cartesian coordinates {X^} in the stress-free state and to coordinates 
{xi} in the deformed state, {i = 1, 2, 3). Uniaxial tension of an incompressible medium is 
described by the formulas 

xi — k(t)X 1 , x 2 = k^{t)X 2l x 3 = k~^(t)X 3 , 

where k — k(t) is the extension ratio. The Finger tensor, F°(t), and the relative Finger 
tensor, F(t,r), read 

F°(t) = A; 2 (t)e 1 e 1 + -^-(e 2 e 2 + e 3 e3), 

k(t) 

nt,r) = (|W) 2 e iei + ||l(e 2 e 2 + e3e 3 ), (30) 

where are base vectors of the frame {Xi}. Substituting expressions (30) into Eq. (29), 
we find the non-zero components (k — 1, 2) of the Cauchy stress tensor 

S = Sidei + S 2 (e 2 e 2 + e 3 e 3 ). 

Excluding the unknown pressure P(t) from the expressions for with the help of the 
boundary condition 

S 2 = 0, 
12 



we find the longitudinal stress 

f f°° r i 

£i(t) = 2p|(l-// L )0Bi(*)+A*L0Li(*)y ex Pry o r(s,w)dsjp(w)ti 
+ [(1 -A*l)^b2(*) +/^l0l 2 (O^ exp (~/ r(s,u;)ds)p(u;)cL 
x [fc 2 (t) + Aj _1 (t)] } [ife 2 (*) - 

+ 2p^ L /;pH^^{0 L1 (t,r) + L2 (t,r)[(M) 2 + ^ 



X 



U(r)/ fc(*) 



r(r, cj) exp — J T(s,uj)ds 



dr. 



(31) 



Equation (31) describes the stress in a specimen for an arbitrary loading program. To 
compare results of numerical simulation with experimental data, we focus on the relax- 
ation test, 

' 1, * < 0, 
A, t > 0, 



k(t) 



(32) 



where A > 1 is a constant. Substitution of Eq. (32) into Eq. (31) implies that 

Si(t) = [<7 B + <r L - <r(t)] (A 2 - A" 1 ), 

where the quantities 

a B = 2p(l-/i L )[0 B1 + 0° 2 (A 2 + A- 1 )\ 
a L = 2p/i L [0° 1 + 0° 2 (A 2 + A- 1 )], 
are independent of time and 

a(t) — <Jl J 1 — exp^— J T(s,cu)ds^j p(u)dw. 

Bearing in mind Eq. (2), we present the dimensionless ratio 

Ei(f) 



in the form 



where 



re 

R(t) = 1 - A 

Jo 



Rit} = E l( 0) 
1 — exp(— r* exp(— uj)t 
A 



p(ui)duj, 



0"L 



(33) 



(34) 



ob + 0" L 

Equations (3) and (33) are determined by 4 adjustable parameters: A is the ratio of the 
relaxing stress to the total stress at the beginning of the test, is the average rate of 
relaxation, whereas Q and £ in Eq. (3) for the function p{uj) describe the distribution 
of relaxation times. These constants are found by fitting experimental data in relaxation 
tests with various elongation ratios, A, for unfilled and particle-reinforced rubbers. 
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9 Experimental 



To assess the time-dependent response of particle-reinforced elastomers, a series of uni- 
axial relaxation tests were carried out at a constant temperature. Dumbbell specimens 
were provided by TARRC laboratories (Hertford, UK) and were used as received without 
mechanical or thermal pre-treatment prior to testing. 

The first compound on the base of natural rubber (EDS-19) contains 1 phr (parts per 
hundred parts of rubber) of high abrasion furnace black (N-330) for UV protection. The 
filler content is too low to affect the mechanical response significantly, and we refer to this 
compound as an unfilled rubber. 

The other compound (EDS-16) contains 45 phr of carbon black (CB) and this com- 
pound is treated as a filled rubber. Before the vulcanization process, a sheet of rubber 
from which the specimens were cut out was milled in a rolling mill to destroy agglomer- 
ates of filler. This procedure leads to a weak anisotropy of the material. The dumbbell 
specimens cut in the direction of milling are referred to as Rl, whereas the samples cut 
in the orthogonal direction are abbreviated as R2. 

A detailed formulation of compounds EDS-16 and EDS-19 is presented in Table 1. 
The mean diameter of CB particles is estimated as 30 nm, and the surface area of filler 
is about 80 m 2 /g [EH . 



Relaxation tests were performed at room temperature at three elongation ratios (A = 
1.2, A = 1.4 and A = 1.8) for samples made of the unfilled rubber, and at four elongation 
ratios (A = 1.2, A = 1.4, A = 1.8 and A = 2.0) for specimens made of the filled rubbers 
Rl and R2. The experiments were carried out by using a testing machine designed at 
the Institute of Physics (Vienna, Austria) and equipped by a video-controlled system. 
To measure the longitudinal strain, two reflection lines were drawn in the central part of 
each specimen before loading (with the distance 7 mm between them). Changes in the 
distance between these lines were measured by a video-extensometer (with the accuracy 
of about 1 %). The tensile force was measured by using a standard loading cell. The 
nominal stress was determined as the ratio of the axial force to the cross-sectional area 
of a specimen in the stress-free state: 4 mm x 2 mm for the unfilled rubber and 4 mm x 
1 mm for the unfilled rubbers. 

To analyze the effects of mechanical pre-loading and thermal recovery, the following 
testing procedure was chosen: 

• Any specimen was loaded with the strain rate e + = 5.0 • 10~ 3 (s _1 ) up to a given 
elongation ratio A, which was preserved constant during the relaxation test (1 hour), 
and unloaded with the strain rate e_ = —5.0 ■ 10~ 3 (s _1 ). 

• The specimen was annealed for 24 hours. 

• The specimen was loaded with the strain rate e + up to the elongation ratio A max = 
3.0 and unloaded with the strain rate e_. 

• The specimen was annealed for 24 hours. 

• The specimen was loaded with the strain rate e + up to the same elongation ratio, 
A, as in the first test. This elongation ratio was preserved constant during the 
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relaxation test (1 hour), and, afterwards, the specimen was unloaded with the strain 
rate e_. 

• The specimen was recovered in an oven for 24 hours at the constant temperature 
T = 100 °C and cooled down by air to room temperature for 4 hours. 

• The specimen was loaded with the strain rate e + up to the same elongation ratio as 
in the first test. This elongation ratio was preserved constant during the relaxation 
test (1 hour), and, afterwards, the specimen was unloaded with the strain rate e_. 

10 Comparison with experimental data 

We begin with matching observations in the first relaxation test on the unfilled rubber with 
the elongation ratio A = 1.2. Because the parameters and Q are mutually dependent 
(an increase in Q results in the growth of r„), we set = 1.0 (s _1 ) in the approximation 
of experimental data. Given a dimensionless ratio A, the quantities Q and E are found 
by the steepest-descent algorithm. The coefficient A is determined by the least-squares 
technique to ensure the minimum discrepancies between experimental data and results of 
numerical simulation. 

We proceed with matching observations in the first relaxation test on the unfilled 
rubber with the elongation ratios A = 1.4 and A = 1.8 by using the parameters Q and 
E found in the approximation of experimental data at the elongation ratio A = 1.2. The 
relaxation rate is determined by the steepest-descent procedure, and the coefficient A 
is found by the least-squares method. 

Figure 1 demonstrates fair agreement between observations and results of numerical 
simulation. It is worth noting that the relaxation curves at different elongation ratios, A, 
noticeably differ from each other, despite some scatter in the experimental data. 

Afterwards, we approximate experimental data in the second relaxation test (after 
pre-loading) on the unfilled specimens with the elongation ratios A = 1.2, A = 1.4 and 
A = 1.8. In the fitting procedure, we employ the values of Q and E which were obtained by 
matching observations in the first relaxation test. The quantity is determined by the 
steepest-descent algorithm, and A is found by the least-squares technique. Experimental 
data and results of numerical analysis are plotted in Figure 2. 

Finally, we fit experimental data in the third relaxation test (after thermal recovery) 
on the unfilled rubber with the elongation ratios A = 1.2, A = 1.4 and A = 1.8. We use 
the same values of Q and E which were employed in matching observations in the first 
and second relaxation tests. The rate of breakage is found by the steepest-descent 
method, and the dimensionless ratio A is determined by the least-squares technique. Fig- 
ure 3 demonstrates good agreement between experimental data and results of numerical 
simulation. 

Using the values of A and found by fitting observations, we calculate the average 
relaxation time, 

n = r-\ 
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and the dimensionless parameter 

A 

r = . 

I — A 

It follows from Eq. (34) that r equals the ratio of the stress in regions with low concen- 
trations of junctions, <7l, to the stress in the bulk medium, <7b, 

r - — 

The relaxation time r* and the coefficient r are depicted in Figures 4 and 5 versus the 
first invariant of the Cauchy deformation tensor, ^(C ). The effect of straining on r* and 
r is described by the linear equations 

r* = r + Ti(h - 3), r = r + ri(/i - 3), (35) 

where the constants (to > 0) and (A; = 0, 1) are determined by the least-squares tech- 
nique. Figures 4 and 5 show that phenomenological equations (35) provide an acceptable 
accuracy of approximation. 

We proceed with matching experimental data for the two grades (Rl and R2) of filled 
rubber by applying the same procedure of fitting. Observations in relaxation tests and 
results of numerical simulation are plotted in Figures 6 to 8 for the samples Rl and in 
Figures 11 to 13 for the specimens R2. The parameters r* and r are depicted versus I\ in 
Figures 9 and 10 for the rubber Rl and in Figures 14 and 15 for the rubber R2 together 
with the results of numerical analysis based on Eq. (35). 

Adjustable parameters for various grades of unfilled and CB filled rubbers are listed 
in Tables 2 to 4, where the abbreviations "V", "P" and "T" are used to denote material 
parameters found in the first (virgin specimens), second (pre-loaded specimens) and third 
(specimens recovered at the elevated temperature) relaxation tests, respectively. 



11 Discussion 

Figures 1 to 3, 6 to 8, and 10 to 12 reveal that the mechanical responses of rubbers (both 
unfilled and particle-reinforced) in the relaxation tests strongly depend on time. The 
relative decrease in the longitudinal stress, 

_ Ei(0)-£i(t ) 

C(to) " — s^o) — ' 

after relaxation for t — 1 hour at the minimum elongation ratio A = 1.2 equals 6 % for 
the unfilled rubber and varies from 13 to 15 % for the CB filled rubbers. The quantity 
C{to) for the filled rubbers (Rl and R2) noticeably exceeds (about by twice) that for the 
unfilled rubber. The fact that the presence of filler enhances relaxation of stresses in 
elastomers may be explained by an increase in the number of regions with low concentra- 
tions of junctions in a particle-reinforced material compared to that in an unfilled medium 
(loading of rubber with the filler particles results in an extra level of inhomogeneity in 
the elastomeric matrix). However, because the dimensionless ratios C(*o) have the same 
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order of magnitude for unfilled and filled rubbers, we suppose that the viscoelastic behav- 
ior of elastomers cannot be described exclusively by means of the mechanically-induced 
destruction of filler clusters and slippage of macromolecules along the surfaces of particles 

10- 



Fair agreement between experimental data in three relaxation tests for the three grades 
of rubber and results of numerical simulation with fixed values of Q and S leads us to the 
conclusion that the relaxation spectrum (which is described in the model by the distri- 
bution function, p(u>), for strands with various energies of breakage, u) is not affected by 
mechanical pre-loading and thermal treatment of specimens. The adjustable parameters, 
Q and E, accept similar values for the two grades of filled rubber (Rl and R2), which 
noticeably differ from those for the unfilled rubber, see Table 2. At the same time, the 
width of the quasi-Gaussian distribution of energies for breakage (which is characterized 
by the dimensionless ratio £ ) is practically independent of the presence of filler (£ belongs 
to the interval between 0.5 and 0.6 for all specimens). The average energy for breakage of 
strands, Q, for the unfilled rubber exceeds by a factor of 2 (approximately) that for the 
filled rubbers. This means that in the domains with low concentrations of junctions where 
breakage and reformation of strands occur, the strength of links between filler particles 
and polymeric chains is essentially weaker than that for chemical and physical cross-links 
and entanglements. 

Figures 2, 6 and 9 demonstrate that the characteristic time of relaxation, r*, increases 
with the longitudinal elongation A. This conclusion contradicts the strain-time superpo- 



sition principle [24, 30], which is conventionally developed within the free- volume theory 
0. According to that concept, mechanical loading of a specimen results in the growth 
of the free volume, which, in turn, accelerates the reformation process for mobile units 
(cooperatively rearranging regions). The free- volume theory correctly describes the time- 
dependent behavior of glassy polymers, where mechanically induced changes in the specific 
volume are substantial, whereas the applicability of this approach to rubbery polymers 
appears to be questionable, because of the incompressibility of rubbers. 

We assume that an increase in the relaxation time, r* (and the corresponding de- 
crease in the relaxation rate, T*) with A may be associated with the mechanically-induced 
crystallization of polymeric chains, see, e.g., [p~0|] . 

In a stress-free state of a rubbery polymer, strands between junctions are composed 
of a number of statistically independent segments freely connected in sequel. For a flex- 
ible chain, the total mechanical energy of these segments is negligible compared to the 
configurational component of the free energy, 

A^ c = -TS C . 

Here S c = /ceOniV) is the configurational entropy per chain, N is the number of con- 
figurations available for a strand with a given end-to-end vector, and the angle brackets 
stand for the mean value. Stretching of a strand causes an increase in the absolute value 
of the end-to-end vector, which results in a substantial decrease in the number of config- 
urations, N. As a consequence, S c decreases and the configurational component of the 
free energy, A\l/ C , increases. At some elongation ratio, A, the free energy for a flexible 
strand exceeds that for a partially crystallized one (where some segments are organized in 
a crystallite structure, whereas the others remain freely jointed), and the crystallization 
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process begins. Because stress-induced crystallization of macromolecules is a kinetic pro- 
cess, it begins from some nuclei of crystallites (segments with restricted mobility). The 
junctions between chains are supposed to play the role of these nuclei, which implies that 
the growth of crystallites occurs from the ends of a strand to its center. 

Crystallization of strands in the neighborhoods of junctions substantially increases 
their strength, which results in a decrease in the rate of breakage, r* (and an increase in 
the relaxation time, r*) with the elongation ratio. 

It is worth noting an important difference between the crystallization processes in 
unfilled and particle-reinforced rubbers. 

In an unfilled elastomer, all junctions in RLCJs are formed by chemical and physical 
cross-links and entanglements. It is assumed that elongation of a strand between these 
junctions results in the growth of crystallites, which practically disappear after unloading 
and subsequent annealing at room temperature. After stretching of the specimens to the 
elongation ratio A = 3.0, the crystallization process during the second relaxation test 
becomes substantially less pronounced, which is concluded from the fact that the slope of 
curve 1 in Figure 4 exceeds that for the curve 2 by a factor of 4. This "training" effect 
(a decrease in the degree of stress-induced crystallization after straining a sample with a 
large amplitude) may be explained by 

1. mechanically- induced alignment of polymeric chains along the direction of stretch- 
ing, see, e.g., 0, 

2. onset of weak links between aligned macromolecules (based on van der Waals forces 
between segments of different chains). 

We suppose that only a part of weak inter-chain links is annihilated when the samples 
are annealed at room temperature, which implies that the mobility of chains in pre-loaded 
specimens (at the second relaxation test) is substantially lower than that in virgin samples. 
As a consequence, the degree of stress-induced crystallization of strands in RLCJs in the 
second relaxation test is noticeably reduced compared to that in the first relaxation test, 
which implies that the relaxation rate, r*, estimated at the second test exceeds that at 
the first test for all elongation ratios. 

Thermal recovery of a (pre-loaded) unfilled rubber at an elevated temperature results 
in (partial) breakage of inter-chain links driven by thermal fluctuations. Thermal treat- 
ment reduces the level of alignment of previously stretched chains in RLCJs and increases 
their ability to be crystallized. This conclusion is confirmed by comparison of curves 2 
and 3 in Figure 4. This figure reveals that for all elongation ratios, A, the relaxation times, 
r* , of the thermally treated specimens exceed those for the untreated samples (which may 
serve as a corroboration for our hypothesis, because the rate of breakage for active strands 
is strongly affected by the degree of their crystallinity). 

For a particle-filled rubber, a part of junctions in a temporary network is formed by 
chemical and physical cross-links, whereas the other part is formed by isolated filler par- 
ticles and small aggregates of filler. We adopt a conventional concept for the rubber-filler 
interaction, see, e.g., ||31|| , according to which filler particles in an elastomeric matrix are 



surrounded by glassy-like cores which are formed at the stages of mixing and vulcaniza- 
tion of rubbery compounds. These hard shells consist of polymeric chains with highly 
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restricted mobility rigidly connected to particles of filler. By analogy with unfilled elas- 
tomers, we suppose that strands bridged to the hard cores are partially crystallized under 
stretching. It is assumed, however, that the hard cores resist the crystallization process, 
which implies that given an elongation ratio, the degree of crystallization in a sample of 
a filled rubber is lower than that in an unfilled one. On the other hand, at unloading 
of a specimen, the glassy-like shells surrounding filler particles prevent annihilation of 
micro-crystallites in active strands linked to these particles, which implies that residual 
crystallites remain in a particle-reinforced elastomer, in contrast to an unfilled rubber, 
where crystallites totally disappear when a specimen returns to its stress-free state. 

The difference between unfilled and filled elastomers may be explained within the 
concept of preferential segregation for polymer chains in the vicinity of filler particles, 
see, e.g., (§. According to this approach, during the composite processing, short chains 
are segregated to interfaces between the elastomeric matrix and particles, whereas the 
concentration of long chains in the host material is increased. This means that short 
chains are mainly bridged to the hard cores. Because the molecular mobility of these 
chains is reduced compared to the long ones, the rate of growth for crystallites in the 
short strands is essentially lower than that in the long strands, which is observed in 
experiments as a very slow increase in the relaxation time, r*, with the elongation ratio, 
A. Comparison of Figure 4 with Figures 9 and 14 implies that the "rate of growth" for 
the relaxation time with the measure of deformation, Ji — 3 [estimated by means of the 
parameter t% in Eq. (35)] for the virgin specimens of the unfilled rubber exceeds that for 
the virgin samples of the filled rubbers by a factor of 30. 

Glassy-like shells surrounding filler particles are assumed to prevent the crystallites 
formed in short strands under stretching from their destruction, which implies that 
stretching a specimen up to A = 3.0 and its subsequent unloading do not affect no- 
ticeably the concentration of crystallites in particle-reinforced rubbers. This is evidenced 
as a pronounced increase in the relaxation times, r*, measured in the second relaxation 
tests on the specimens of filled rubbers compared to that in the relaxation tests on the 
virgin samples, see Figures 9 and 14. 

Thermal recovery of a filled rubber at an elevated temperature causes destruction of 
crystallites arisen under loading. This process is supposed to be driven by (i) thermal 
fluctuations which substantially accelerate breakage and reformation of chains linked to 
the filler particles and (ii) thermal annihilation of glassy-like cores around particles. As 
a result, the degree of crystallinity in reinforced specimens is reduced after recovery, and 
the relaxation rate, r*, is increased. The degree of crystallinity estimated in the third 
relaxation test exceeds, however, that for virgin specimens (see curves 1 and 3 in Figures 
9 and 14), which means that some residual crystallites remain even after 24 hours of 
recovery at T = 100 °C. 

Table 3 evidences that the parameter tq for the filled rubber R2 (samples cut out 
in the direction orthogonal to the direction of milling) noticeably exceeds that for the 
filled rubber Rl (specimens cut in the direction of milling) in all three relaxation tests. 
This may be explained by alignment of chains in RLCJs during the milling procedure. 
Polymeric chains with the end-to-end vectors directed along the direction of loading in 
the relaxation tests are strongly aligned during milling, which substantially reduces their 
mobility (because of the presence of inter-chain links formed at alignment) and, as a 
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consequence, their ability to be crystallized. Because these chains bear the main load 
at uniaxial extension of specimens, the average relaxation rate for specimens Rl remains 
rather high. On the contrary, alignment of polymeric chains with the end-to-end vectors 
directed perpendicular to the direction of stretching for specimens R2 unsubstantially 
affects the relaxation rate (these chains are weakly strained in tensile relaxation tests, and 
their mechanical response is negligible), whereas the chains directed along the stretching 
direction are practically not aligned, and the degree of their crystallization is similar to 
that for the chains in a specimen not subjected to milling. 

Figures 5, 10 and 15 demonstrate that for all grades of rubber the dimensionless ratio 
r increases with the longitudinal strain in the first and third relaxation tests and decreases 
in the second test. This means that in the first and third tests, the stress, <tl, in RLCJs 
grows more rapidly (with the elongation ratio A) than the stress, ctb, in the bulk material, 
while in the second test, the stress in the bulk medium increases with A substantially 
more rapidly than the stress, <Tl, in domains with low concentrations of junctions. This 
qualitative difference between the mechanical response in the first and third tests, on 
the one hand, and in the second test, on the other, is observed for both unfilled and 
particle-reinforced rubbers. It may be associated with the stress-induced alignment of 
chains. 

Experimental data on virgin samples show that uniaxial tension of rubber causes a 
more pronounced growth of the stress, <tl, in RLCJs than the growth of the stress, <tb, 
in the bulk material. This difference may be explained by the mechanically-induced 
crystallization of active strands in the regions with low concentrations of junctions (we 
presume that stress-induced crystallization of chains either does not occur in the bulk 
medium or takes place with a substantially lower rate than in RLCJs). 

After pre-loading to A = 3.0 and subsequent annealing at room temperature, the 
stress in the bulk material grows with the elongation ratio, A, more rapidly than the 
stress in RLCJs. This observation is associated with a higher level of alignment of chains 
in the bulk material (which causes its strain-hardening) than that in domains with low 
concentrations of junctions. We assume that breakage and reformation of active strands 
during annealing for 24 hours noticeably destroy alignment of chains in RLCJs driven 
by the mechanical pre-loading, whereas alignment of chains in the bulk material (where 
reformation of strands does not take place) remains at the same level as immediately after 
pre-loading. 

Finally, after thermal recovery at an elevated temperature, the stress in RLCJs in- 
creases with the elongation ratio more rapidly that the stress in the bulk material, which 
is explained (as in the case of virgin specimens) mechanically-induced crystallization of 
strands in inclusions with low concentrations of junctions. 

It is worth noting that the dimensionless ratio r for the unfilled rubber accepts the 
same values for the first and third relaxation tests (because crystallites disappear in an 
unfilled elastomer after unloading and thermal recovery), while the ratio r for the third 
test on filled rubbers exceeds that for the first test (which confirms our assumption that 
some residual crystallites remain in a particle- reinforced elastomer after unloading). 

Comparison of Figure 5, on the one hand, and Figures 10 and 15, on the other, reveals 
that the dimensionless parameter r for the unfilled rubber exceeds (about by twice) that 
for the filled rubbers. This fact may be explained by a higher level of inhomogeneity of 
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particle-reinforced elastomers compared to the unfilled ones, which implies that RLCJs 
in filled rubbers are noticeably "looser" than those in unfilled rubbers with the same 
polymeric matrix. 

Fitting of experimental data implies practically the same values of the dimensionless 
ratio r for specimens Rl and R2, which means that the effect of milling on this parameter 
is negligible. 

12 Concluding remarks 

Constitutive equations have been derived for the isothermal viscoelastic behavior of filled 
elastomers at finite strains. With reference to the concept of temporary networks, a filled 
elastomer is treated as a network of long chains connected to junctions (chemical and 
physical cross-links, entanglements and filler particles). It is assumed that a particle- 
reinforced rubber is a strongly non-homogeneous medium (with the characteristic length 
of inhomogeneity of about 1 /im) which consists of an elastomeric matrix with a normal 
concentration of cross-links, wherein inclusions with low concentrations of junctions are 
distributed. In the domains with low concentrations of cross-links, the average number of 
junctions between polymeric chains (per unit mass) is less than that in the bulk material 
because of local inhomogeneities in the distribution of a cross-linker at mixing. 

It is assumed that intermolecular forces in the bulk material prevent slippage of chains 
from junctions (which are treated as permanent), whereas in RLCJs active strands (whose 
ends are linked to separate junctions) can slip from the temporary junctions as these 
strands are thermally agitated. Breakage of an active strand (its transition into a dangling 
one) occurs at a random instant when one of its ends detaches from a junction. An active 
strand is reformed at a random time when the free end of a dangling strand captures some 
junction in its vicinity. The stress in a dangling strand totally relax before this strand 
merges with the network, which implies that the stress-free state of an active strand 
coincides with the deformed state of the network at the instant of their merging. Unlike 
conventional theories of transient networks, we postulate that different strands in RLCJs 
have different energies for breakage and describe these energies by the quasi- Gaussian 
distribution function, p(u), with two adjustable parameters, fl and E. 

Constitutive equations for a filled rubber are developed by using the laws of thermo- 
dynamics. These equations are applied to determine stresses at uniaxial extension of a 
specimen. The time-dependent response of a filled rubber in the standard relaxation test 
is characterized by four adjustable parameters which are found by fitting observations. To 
determine these constants, a series of relaxation tests have been carried out on specimens 
of unfilled and CB filled natural rubber (virgin samples, samples pre-loaded up to the 
elongation ratio A = 3.0, and specimens pre-loaded and recovered at an elevated tem- 
perature). Fair agreement is demonstrated between results of numerical simulation and 
observations in relaxation tests (with the duration of t — 1 hour) at various elongation 
ratios, A, in the region of moderate finite deformations, A G (1,2). 

The following conclusions are drawn: 

1. The unfilled and particle reinforced rubbers demonstrate a pronounced viscoelastic 
response. The relative decrease in stresses, ((t ), for the filled rubbers exceeds that 
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for the unfilled rubber by twice. This difference may be ascribed to an increase in 
the number of regions with low concentrations of junctions caused by the presence 
of filler. 

2. The average energy for breakage of strands, Q, decreases with an increase in the filler 
content, whereas the width of the quasi- Gaussian distribution of energies, £ = E/f2, 
is independent of the content of carbon black. The distribution function for the 
energies of breakage, p(u), is not affected by strains, mechanical pre-loading (up to 
the elongation ratio A = 3.0) and thermal recovery. 

3. The characteristic time of relaxation, r*, increases with the elongation ratio, A. 
The growth of r* with A is explained by the mechanically-induced crystallization of 
natural rubber in the regions with low concentrations of junctions. 

4. A qualitative difference is revealed in the influence of mechanical pre-loading and 
thermal recovery on the relaxation time, r*, for unfilled and particle-reinforced elas- 
tomers. The characteristic time of relaxation for the unfilled rubber is substantially 
reduced after pre-loading, and it weakly increases after thermal recovery at an el- 
evated temperature. On the contrary, mechanical pre-loading of the filled rubber 
causes a noticeable growth of its relaxation time, which decreases after thermal re- 
covery. This difference is ascribed to the influence of glassy-like cores around filler 
particles, which prevent annihilation of crystallites in active strands bridged to the 
particles of filler. 
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Table 1: Chemical composition of rubbers (phr by weight) 



Formulation 


EDS-16 


EDS-19 


Natural rubber 


100 


100 


Zinc oxide 


5 


5 


Stearic acid 


2 


2 


Carbon black (N330) 


45 


1 


Process oil 


4.5 





Antidegradant (HPPD) 


3 


3 


Antiozonant wax 


2 


2 


Accelerator (CBS) 


0.6 


0.6 


Sulphur 


2.5 


2.5 



Table 2: Adjustable parameters Q and S 



Rubber 


Q S £ 


Unfilled 

Rl 

R2 


13.80 7.00 0.51 
5.30 2.80 0.53 
5.13 3.18 0.62 



Table 3: Adjustable parameters r and t\ 



Rubber 


To 


Tl 




V P T 


V P T 


Unfilled 

Rl 

R2 


0.0 1.1866 1.9068 
1.0114 1.7418 1.2543 
1.3129 2.8106 1.9364 


10.4615 2.5193 3.1182 
0.1217 0.3390 0.1960 
0.3073 -0.0226 0.1598 


Table 4: Adjustable parameters r and r\ 


Rubber 


To 


n 




V P T 


V P T 


Unfilled 

Rl 

R2 


0.3754 0.4674 0.3754 
0.2008 0.2415 0.2647 
0.2132 0.2655 0.2389 


0.1964 -0.0354 0.1964 
0.0364 -0.0286 0.0310 
0.0416 -0.0453 0.0495 
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Figure 1: The dimensionless ratio R versus time t (s) for unfilled rubber (virgin samples) 
in a relaxation test with the elongation ratio A. Symbols: experimental data. Solid lines: 
results of numerical simulation. Curve 1: A = 1.2; curve 2: A = 1.4; curve 3: A = 1.8 

Figure 2: The dimensionless ratio R versus time t (s) for unfilled rubber (after mechanical 
pre-loading) in a relaxation test with the elongation ratio A. Symbols: experimental data. 
Solid lines: results of numerical simulation. Curve 1: A = 1.2; curve 2: A = 1.4; curve 3: 
A = 1.8 

Figure 3: The dimensionless ratio R versus time t (s) for unfilled rubber (after pre- 
loading and thermal recovery) in a relaxation test with the elongation ratio A. Symbols: 
experimental data. Solid lines: results of numerical simulation. Curve 1: A = 1.2; curve 
2: A = 1.4; curve 3: A = 1.8 

Figure 4: The characteristic time of relaxation r* (s) versus the first invariant I± of the 
Cauchy deformation tensor for unfilled rubber in relaxation tests. Symbols: treatment 
of observations. Unfilled circles: a virgin specimen; filled circles: the specimen after 
mechanical pre-loading; asterisks: the specimen after pre-loading and thermal recovery. 
Solid lines: approximation of the experimental data by Eq. (35). Curve 1: r = 0.0, 
ri = 10.4615; curve 2: r = 1.1866, n = 2.5193; curve 3: r = 1.9068, n = 3.1182 

Figure 5: The dimensionless parameter r versus the first invariant of the Cauchy de- 
formation tensor l\ for unfilled rubber in relaxation tests. Symbols: treatment of obser- 
vations. Unfilled circles: a virgin specimen; filled circles: the specimen after mechanical 
pre-loading; asterisks: the specimen after pre-loading and thermal recovery. Solid lines: 
approximation of the experimental data by Eq. (35). Curves 1 and 3: r = 0.3754, 
r l = 0.1964; curve 2: r = 0.4674, r 1 = -0.0354 

Figure 6: The dimensionless ratio R versus time t (s) for CB filled rubber (Rl) in 
a relaxation test with the elongation ratio A. Circles: experimental data. Solid lines: 
results of numerical simulation. Curve 1: A = 1.2; curve 2: A = 1.4; curve 3: A = 1.8; 
curve 4: A = 2.0 

Figure 7: The dimensionless ratio R versus time t (s) for CB filled rubber (Rl) in 
a relaxation test (after mechanical pre-loading) with the elongation ratio A. Circles: 
experimental data. Solid lines: results of numerical simulation. Curve 1: A = 1.2; curve 
2: A = 1.4; curve 3: A = 1.8; curve 4: A = 2.0 

Figure 8: The dimensionless ratio R versus time t (s) for CB filled rubber (Rl) in 
a relaxation test (after pre-loading and thermal recovery) with the elongation ratio A. 
Circles: experimental data. Solid lines: results of numerical simulation. Curve 1: A= 1.2; 
curve 2: A = 1.4; curve 3: A = 1.8; curve 4: A = 2.0 

Figure 9: The characteristic time of relaxation r* (s) versus the first invariant Ii of the 
Cauchy deformation tensor for CB filled rubber (Rl) in relaxation tests. Symbols: treat- 
ment of observations. Unfilled circles: a virgin specimen; filled circles: the specimen after 
mechanical pre-loading; asterisks: the specimen after pre-loading and thermal recovery. 
Solid lines: approximation of the experimental data by Eq. (35). Curve 1: t = 1.0114, 
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ri = 0.1217; curve 2: r = 1.7418, n = 0.3390; curve 3: r = 1.2543, n = 0.1960 

Figure 10: The dimensionless parameter r versus the first invariant of the Cauchy de- 
formation tensor I± for CB filled rubber (Rl) in relaxation tests at room temperature. 
Symbols: treatment of observations. Unfilled circles: a virgin specimen; filled circles: 
the specimen after mechanical pre-loading; asterisks: the specimen after pre-loading and 
thermal recovery. Solid lines: approximation of experimental data by Eq. (35). Curve 
1: r = 0.2008, n = 0.0364; curve 2: r = 0.2415, n = -0.0286; curve 3: r = 0.2647, 

n = 0.0310 

Figure 11: The dimensionless ratio R versus time t (s) for CB filled rubber (R2) in a 
relaxation test with the elongation ratio A. Circles: experimental data. Solid lines: results 
of numerical simulation. Curve 1: A = 1.2; curve 2: A = 1.4; curve 3: A = 1.8; curve 4: 
A = 2.0 

Figure 12: The dimensionless ratio R versus time t (s) for CB filled rubber (R2) in 
a relaxation test (after mechanical pre-loading) with the elongation ratio A. Circles: 
experimental data. Solid lines: results of numerical simulation. Curve 1: A = 1.2; curve 
2: A = 1.4; curve 3: A = 1.8; curve 4: A = 2.0 

Figure 13: The dimensionless ratio R versus time t (s) for CB filled rubber (R2) in 
a relaxation test (after pre-loading and thermal recovery) with the elongation ratio A. 
Circles: experimental data. Solid lines: results of numerical simulation. Curve 1: A = 1.2; 
curve 2: A = 1.4; curve 3: A = 1.8; curve 4: A = 2.0 

Figure 14: The characteristic time of relaxation r* (s) versus the first invariant I x of the 
Cauchy deformation tensor for CB filled rubber (R2) in relaxation tests. Symbols: treat- 
ment of observations. Unfilled circles: a virgin specimen; filled circles: the specimen after 
mechanical pre-loading; asterisks: the specimen after pre-loading and thermal recovery. 
Solid lines: approximation of the experimental data by Eq. (35). Curve 1: r = 1.3129, 
ri = 0.3073; curve 2: r = 2.8106, n = -0.0226; curve 3: r = 1.9364, n = 0.1598 

Figure 15: The dimensionless parameter r versus the first invariant of the Cauchy de- 
formation tensor Ii for CB filled rubber (R2) in relaxation tests at room temperature. 
Symbols: treatment of observations. Unfilled circles: a virgin specimen; filled circles: 
the specimen after mechanical pre-loading; asterisks: the specimen after pre-loading and 
thermal recovery. Solid lines: approximation of the experimental data by Eq. (35). Curve 
1: r = 0.2132, n = 0.0416; curve 2: r = 0.2655, n = -0.0453; curve 3: r = 0.2389, 
r 1 = 0.0495 
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